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Mixing of chemical species in reacting flows 


Chemically Reacting Laminar Flows 

Numerical (or semi-analytical) methods 

O.D.E. solvers 

1-D Numerical methods 

Numerical uncertainty is minimal 

Chemical kinetics 

Chemical kinetics (finite rate chemical kinetics & Reduced finite rate chemical 
kinetics) for certain fuels (hydrogen and some of hydro-carbons) have been 
tested and validated repeatedly. 

Numerical and theoretical analyses yield accurate results. 


Chemically Reacting Turbulent Flows 


Numerical methods 

Boundary layer or Navier-Stokes equations solvers 
(Except for semi-analytical analyses of PSR cases) 

Uncertainty is caused partly by numerical methods. 

Chemical kinetics 

Finite rate chemical kinetics & Reduced finite rate chemical are used. 

Uncertainty is caused partly by over-simplified chemical kinetics. 

Turbulence Equations 

Numerical analysis yields not so accurate results. Unceratinlty caused by turbulence 
equations is by far greater than that caused by chemical kinetics.* 

The probability density functions are used to improve the predictive capability of 
turbulent mixing of chemical spedes. However, some pdf methods yield physically 
incorrect numerical results. 

It is certainly important to better understand the turbulent mixing of chemical 
spedes in reacting flows. 

♦Westbrook and Dryer, Prog . Energy Combustion Science, vol. 10, pp. 1-57, 1984. 



Turbulent Mixing of Chemical Species in Reacting Flows 


1. Large eddy mixing 

2. Turbulent Mixing 

3. Molecular Mixing 

The use of accurate chemical kinetics alone can not yield accurate results for 
numerical analysis of turbulent reacting flows unless turbulent mixing is resolved 
correctly. 


Large Eddy Mixing 

Caused by separated flows, recirculating flows, and organized structures. 

The extent of calculated recirculation zone or the organized structures depends 
on the accuracy of both the numerical method and the turbulence equations since 
they are nonlinearly coupled with each other. 



W here R = W(k, £t, • • • ), and Pr^j^) 2 * (| + gjf} 

Therefore, the capability to correctly resolve the large eddy mixing depends on 
the accuracy of the numerical method and the turbulence equations. 


Numerical investigations carried out during last decades show that k-£, ARSM, 
and RSM do not yield accurate results for complex turbulent flows. 

On the other hand, the multiple-time-scale turbulence equations* yield highly 
improved numerical results for various complex turbulent flows. 


*Kim and Benson, Int. J. Heat Mass Transfer, vol.35, pp. 2357-2365, 1992. 


Turbulent Mixing 


Caused mostly by energy-containing eddies and partly by fine-scale eddies. 

The capability to resolve the turbulent mixing of chemical species depends on the 
capability of the turbulence equations to correctly describe the turbulent 
transport of scalar variables (i.e., heat transfer, turbulent kinetic energy, 
concentrations, and convection-diffusion of species, etc) 

Need to be able to describe the chemical reaction-turbulence interaction (i.e., 
turbulent mixing is enhanced and shear layer thickness is widened by chemical 
reaction)* 

The multiple-time-scale turbulence equations can resolve the cascade of 
turbulent kinetic and the nonequilibrium turbulence phenomena. 

(Present numerical results show that the M-S turbulence equations can resolve 
the chemical reaction-turbulence interaction.) 


* C.T. Chang et al., "Comparison of reacting and non-reacting shear layers at a high 
subsonic mach number," AIAA paper 93-2381, 1993. 


Molecular Mixing 
Caused by molecular diffusivity. 

Theoretical and numerical analyses of chemically reacting laminar flows yield 
accurate results. 

Various molecular diffusion equations can accurately describe the molecular 
diffusion of species. 

The Lennard-Jones 12-6 potential law is used in the present study. 

Various numerical methods that yield accurate results for laminar flows can 
resolve the molecular mixing. 


Multiple-time-scale turbulence equations 
Energy containing eddies 
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Turbulent Eddy Viscosity 

Pt= pCpk 2 /e p where = c^ f et/e p 

Remark: Single-time-scale turbulence models can not resolve nonequilibrium 
turbulence phenomena. 


Species conservation equations for reacting flows 
Chemical species concentration equation 

^(pYi) + V- {pYi(v + Vi)}= wj 

where the diffusion velocity, Vj, is approximated using the Fick's law given as 
YiV i = -(D iJ + D i(t )VY i 

The production rate of the i-th species, wj, is given as 



Chemical reactions for the combustion of H 2 in a vitiated supersonic airstream 
are described using 9 chemical species (H 2 , O 2 , H 2 O, OH, O, H, HO 2 , H 2 O 2 , and 
N 2 ) and 24 pairs of reaction-steps (Burks and Oran, 1981; Kumar, 1989). 

A fast chemistry can not be used to describe the fine details of chemically 
reacting flows. 

A reduced chemical kinetics can not be used confidently due to the uncertainty 
contained in the reaction mechanisms. 

The use of a detailed finite rate chemistry may make it difficult to obtain a fully 
converged solution due to the coupling between the large number of flow, 
turbulence, and chemical equations. The numerical method needs to be 
strongly convergent. Accuracy also depends on the capability of turbulence 
equations used. 


Numerical Method 

The numerical method is a finite volume method that incorporates a 
pressure-staggered mesh and an incremental pressure equation for the 
conservation of mass. 

Predictor Step: Solve momentum equation. 

(pC 1+ A*)u**=I A &*+ Sj- + pC 2 uf >- pC 3 un- 2 (1) 
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Corrector Step: Correct the velocity field to be divergence free 
Incremental pressure equation 
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Incremental velocity equation 


Velocity and pressure corrections 


U*** = u*‘+u'i 

(4) 

p" = p* + p' 

(5) 


Solve eqs. (1-5) iteratively until all flow variables are converged. 



Combustion of H 2 in vitiated supersonic airstream (Burrows and Kurkov, 1973) 
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Calculated OH contour. 


The numerical result obtained using the M-S turbulence equations indicates 
that the combustion occurs at approximately middle of the channel. The 
calculated flame location is in excellent agreement with that observed in the 
experiment. 




Temperature contours obtained using the M-S turbulence equations 


It can be seen in the temperature contours that the calculated flame location 
obtained using the M-S turbulence equations are in correct agreement with the 
measured data. The contours also indicate that the temperature increases 
within in a short distance. The trend is in correct agreement with experimental 
observations that temperature increase occurs within a finite flame thickness. 



Temperature contours obtained using the pdf* 


The pdf fails to predict the correct flame location (i.e., ignition delay). The 
slowly increasing temperature field indicates that the pdf may not be able to 
predict a correct flame front. The numerical results obtained using the pdf are 
not in correct agreement with the physics of combustion. 


* A. T. Hsu, Y.-L.P. Tsai, and M. S. Raju, "A PDF approach for compressible 
turbulent reacting flows," AIAA Paper 93-0087, 1993. 




Conclusions and Discussion 


The calculated species concentration profiles are in as good agreement with 
the measured data as those obtained using the pdf. 

The flame location (ignition delay) obtained using the multiple-time-scale 
turb ulence equations is in excellent agreement with the experimentally 
observed onset of ultraviolet radiation. 

Cascade of the turbulence filed is influenced by the extra strains caused by 
chemical reaction. 

Both the numerical results and the measured data exhibit enhanced mixing 
of the hydrogen and vitiated airstream for the reacting case. 

The M-S turbulence equations can resolve the chemical reaction-turbulence 
interaction. 

The pdf produces slowly increasing temperature field and it fails to predict 
the ignition delay. Thus the numerical results obtained using the pdf are not in 
correct agreement with physics of combustion. 


Part B. Unsteady Transitional Flows over Forced Oscillatory Surfaces* 
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* S.-W. Kim, K.B.M.Q. Zaman, and J. Panda, "Calculation of unsteady transitional 
flow over oscillating airfoil" in Separated Flows, eds. J.C. Dutton and L.P. Purtell, 
Proceeding of ASME Fluid Engineering Conference, Washington D.C., June 20-24, 
1993. 



Nomenclature 



Lagrangian-Eulerian coordinates 

x: fixed reference coordinates 
moving coordinates 
J = I3£j/3ajl 


Unsteady transitional flow equations with moving boundaries 


Conservation of mass equation 

di^\ T 3 


a) 


Conservation of linear momentum equation 
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Numerical Method 


The unsteady transitional flow equations are solved using the same finite 
volume method. Time-integration is made using an iterative-time-advancing 
scheme. 

Comparison of Unsteady Flow Solution Techniques* 

Iterative Time-Advancing Scheme (ITA). 

Simplified Marker and Cell (SMAC). 

Pressure-Implicit Splitting of Operators (PISO). 

The ITA that can best resolve the nonlinearity of the Navier-Stokes equations, 
and yields the most accurate results. 

The SMAC is the most efficient computationally and yields accurate numerical 
results for laminar flows. 

The PISO is the most unstable numericaly and yields less accurate results. 


* Kim & Benson, Computers and Fluids, vol. 21, pp. 435-454, 1992 



(a) a = 15°' t/T = 0.0 
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Oscillating airfoil and moving mesh 
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Ensenble-averaged velocity profiles at x/ c=1.0 
The deteriorated comparison at a =19°d 

(i) The hot wire can not accurately measure the velocity components when the 
flow is misaligned more than approximately 30° from the hot wire axis. 

(ii) The interaction between the DSV and the TEV occurs in a relatively coarse 
mesh region and the numerical method yields somewhat deteriorated results 








Turbulent viscosity contour ( A|i/p t = 10) 



CONCLUSIONS AND DISCUSSION 

Numerical method successfully predicts the Dynamic Stall Vortex and the 
Trailing Edge Vortex. 

The calculated ensemble-averaged velocity profiles are in good agreement 
with the measured data. 

Both the numerical results and the measured data show that the transition 
from laminar to turbulent state and relaminarization occur widely in space and 
in time. 

The good comparison between the numerical results and the experimental 
data are attributed to the capability of 

(i) the ITA that can best resolve the nonlinearity of the Navier-Stokes 
equations, 

(ii) the new pressure correction algorithm that can strongly enforce the 
conservation of mass, and 

(iii) the Multiple-time-scale turbulence equations that can resolve the 
transitional nonequilihrium turbulence field. 


